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Abstract 

By coupling the asymmetric three-terminal mesoscopic dielectric system with a temperature probe, at 
low temperature, the ballistic heat flux flow through the other two asymmetric terminals in the nonlinear 
response regime is studied based on the Landauer formulation of transport theory. The thermal rectification 
is attained at the quantum regime. It is a purely quantum effect and is determined by the dependence 
of the ratio trc(co)/trl(u>) on a>, the phonon's frequency. Where trc{(jS) and trl(u>) are respectively the 
transmission coefficients from two asymmetric terminals to the temperature probe, which are determined by 
the inelastic scattering of ballistic phonons in the temperature probe. Our results are confirmed by extensive 
numerical simulations. 
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At low temperature, the thermal conductance in mesoscopic phonon systems is dominated by 
the transmission of ballistic phonons. Using the Landauer formulation of transport theory, Rego 
and KirczenowyJ] predicted the existence of the quantum of ballistic thermal conductance. Fur- 
ther investigations indicated that the quantum of ballistic thermal conductance is universal and it 
is independent of the statistics obeyed by the carriers [2- 4]. The quantum of ballistic thermal con- 
ductance has already been verified by experiment[5] and the universality is verified later[6]. Since 
then, the quantum of ballistic thermal conductance has attracted much attention [3- 11]. Besides 



the unique property of the quantized thermal conductance in the mososcopic systems, the other 
motivation for these investigations is to seek ways to control the heat transport at low temperature. 

To control the heat transport, the first nanoscale model for thermal rectifier was proposed in 
ID nonlinear chains in 2002 111211 . After that, many theoretical works were performed to study the 



thermal rectification|l3 



21Q . and the first experimental work to demonstrate the thermal rectifica- 



tion was also reported|22Q, although only the 2 to 7% rectification efficiency is yielded. Recently, 
based on the ID nonlinear chains, the models of thermal logic gatesj23] and thermal memory[24] 
were proposed to demonstrate that phonons can be used to carry information and processed ac- 
cordingly. Although considerable achievements in thermal rectification based on the nonlinear 
crystal are made, the thermal rectification based on the harmonic system is still unavailable. This 
may due to the fact that there is no thermal rectification in the heat transport in a two-terminal har- 
monic system at low temperature QJ/ZD- This absence of thermal rectification in harmonic system 



26Q and 



is also found in classical models: the harmonic chains with self-consistent reservoirs [25, 
harmonic /V-terminal j unction Il27ll . However, at the quantum regime, for an asymmetric harmonic 
system which couples with a temperature probe connecting directly to a self-consistent reservoir, 
the temperatures of the self-consistent reservoir for forward and reversed bias are not symmetric 
with respect to the average temperature[28]. Therefore, one can expect that the thermal rectifica- 
tion can be attained in this asymmetric harmonic system at the quantum regime. 

In this work, we theoretically demonstrate the possibility to attain the ballistic thermal rectifi- 
cation by coupling the asymmetric mesoscopic three-terminal junctions with a temperature probe. 
This possibility is confirmed by extensive numerical simulations. Just same as the voltage probe 



in the mesoscopic electronic systems [2 



which has been used to study the electrical rectifica- 

n 



tion in three-terminal electrical ballistic junctions |31|], the temperature probe acts as a dephasing 
probe which introduces inelastic scattering of phonons in three-terminal junction. This inelastic 
scattering brings nonlinearity in the asymmetric three-terminal junctions and makes the thermal 



rectification possible. 

The asymmetric three-terminal junction is sketched in the inset of Fig. [TJ Through terminals L 
and C, the junction is coupled to two thermal reservoirs at temperatures T L and T c . The steady- 
state heat flux Q = Qc = —Ql is passed through the junction via terminal C and L. A third terminal 
R is connected to another thermal reservoir at temperature T R . This third terminal is a temperature 
probe, it means T R is adjusted in such a way that no net heat flux pass through terminal R. The 
ballistic regime is considered in this work. The heat flux <2, from terminal i (i = L, R, C) flowing 



into the midsection J can be expressed as [7, 



32] 



ZC 00 dco 
[n(T h to) - n(Tj, co)\hco Tji(co)—, (1) 

where n(Tj,co) = [exp(hco/k B Ti) - 1]~ ! is the Bose-Einstein distribution function of the phonons 
in the zth reservoir, T t is the equilibrium temperature of thermal reservoir i, t ;; (cl>) = Yim _ 
0Jim)T ji, m (a>) = Dm,n _ u j n )6(to - w, ot )t ji,nm(co) is the total transmission coefficient and r ;i -,„ m (ct>) 
is the transmission coefficient from mode m of terminal i at frequency co across all the interface into 
the mode n of terminal j, oj im is the cutoff frequency of mode m in terminal i. By the time-reversal 
symmetry, we have t^&j) = Tij{co). The influence of the temperature probe R on the steady-state 
heat flux Q is determined by the transmission coefficients trc(co) and trl(co). 

Firstly, the reservoir L is set at temperature T L = T/ 10t , which is higher than the temperature 
Tc = T co id of reservoir C. The heat flux flows through the system is <2_ = Qc- Evidently, the heat 
flux flows into C reservoir hence Q < 0. By reversing the temperature bias, i.e., let T L = T co u 
and T c = That, the heat flux equals to Q+ = Qc > 0. 

To study whether the ballistic thermal rectification will exhibit, <2_ is added to Q + . 



[n(T C0 



+n(T hot , 0) - n(T' R , io)]T RC {co)fico (2) 



AQ = + (J = I [n(T coId , co) - n(T R , co) 

dco 
lit 



where T R and T' R are determined by Q R = 0. Obviously, T R and T' R are determined by the transmis- 
sion coefficients t R c(co) and t R l(co). Therefore, whether the rectification can be attained depends 
critically on the transmission coefficients r RC (co) and t rl (co). In the linear response regime, by 



Qi = ^^(7,-^)10, 



using (ji = 2^j±i trji^ i - 1 j ) IL/J, 1220, where G ; , is the two-terminal thermal conductance from reser- 
voir i to j, one can easily find that n(T cokl , co) - n{T R , co) + n(T hot , co) - n{T' R , co) = for any co. Thus 
AQ = 0, which means there isn't rectification in the linear response regime. This corresponds 



to the one obtained in 03311 for the thermoelectric transport in a chain of quantum dots with self- 
consistent reservoirs. On the other hand, if the ratio trc(co)/trl(co) is a constant for any frequency 
co, by using Q R = 0, one can easily find AQ = and the rectification is absent even in the nonlinear 
response regime. Therefore, even in the nonlinear regime, the thermal rectification is absent in the 
symmetric three-terminal junctions with t rc (co) = t rl (co). 

In a more general situation where the ratio Tr C (co)It rl (co) for any frequency is not a constant, 
the fact that trc(co)/trl(co) varies with co means there are different asymmetries for transmission 
of phonons in different frequencies. One can certainly expect that AQ is not zero in the nonlinear 
response regime and it is determined by the dependence of the ratio t rc (co)/tr L (co) on co. We 
start with letting T = (T hot + r coW )/2 and AT = (T hot - T co i d )/2, where T can be regarded as 
the working temperature according to the following text. When AT is a small value, by using the 
Taylor expansion of the Bose-Einstein distribution function n(Ti, co) as well as the approximate 



results of T R and T' R in ref. Il28h for the asymmetric three-terminal junctions, one can find[34] 



_ ^ 02 ^2(Trc)Fi(Tr L ) - F 2 (Tr L )F 1 (Tr C ) ^ t ^2 
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+<9[(Ar) 4 ] = a(AT) 2 + 0[(AT) 4 ]. (3) 
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By a simple algebra, F 2 {trc)Fi(t rl ) - F 2 (tr L )Fi(trc) can be reexpressed as 



where 
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T(C0 u a) 2 ) = -T(C0 2 ,C0y) = Tr C (C0 y )t rl (C02) ~ Tr C (C02)Tr L (C0 x ). 



(6) 



(7) 
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This is the main result of this work. In Eq. ©, d 2 n/dT 2 = [(1 + 2n)hco I (k B T 2 ) - 2/T]dn/dT 
is used. Obviously, N(co\,co2, Tq) is T dependent and is originating from the quantum statistics 
of phonons. T(co\, to 2 ) is originating from the inelastic scattering of phonons in the temperature 
probe. The first term of T(to\, CO2) in Eq. © can be understood as that the ballistic phonon with 
frequency tO\ transport from reservoir C into temperature probe R with transmission coefficient 
Tr C (co\), after being inelastically scattered into the different frequency co 2 , the ballistic phonon 
transport from the temperature probe into reservoir L with transmission coefficient t rl (co 2 ). The 
second term presents the similar process but in the inverse direction. Thus the ballistic thermal rec- 
tification is a purely quantum effect and is determined by the inelastic scattering in the temperature 



probe. To study the sign of N{co\,co 2 , 7 ), d{x[l + 2n(x)]}/dx is studied, where x = ~hcolk B T . 

d{x[\ + 2n(x)]} e 2x - 2xe x - 1 



2 



dx (e x - 1) 

for all x > 0. Hence N(coi,co2,T ) > for all co\ > a>2 > 0. Thus, the inelastic scattering 
factor, T(co\, ojz), determines the sign of a due to < 1 < 1. By reexpressed T{co\, C02) as 
[t sc (^»i)/trl(^i) - t rc (co2)It rl (<x>2)\t rl (cox)t rl (<jJ2) when r«i(w) * for all co, one can easily 
find that the sign of a is determined by the dependence of the ratio t rc {co)It rl {co) on co. AQ > 
when t rc {co)It rl {co) increases with co. 

To obtain the exact AQ, we carry out numerical calculations for Eq. © with the scattering 



matrix method used in [28]. In the calculation, the geometrical parameters of the system are 
chosen as W L = W R = W c = Dj = 10 nm while W T can be varied. We limit the temperature T hot 
of the thermal reservoir at the higher temperature is lower than T ph = ~hnv/W L k B ~ 7.61 K (v is the 
sound velocity) to ensure that the phonon relaxation can be neglected [7] and the heat conduction 
is determined by the ballistic transmission of the acoustic phonons. 

Fig. Q] shows the results of AQ vs 2A7 for different 7 with Wj = 21 nm. Firstly, at small 
but finite A7, AQ shows a quadratic dependence on A7, in agreement with Eq. ©. Secondly, 
corresponds to the factor N(u>i,a>2,To), AQ decreases with increasing 7 for a system with a 
certain W T . Thirdly, when the temperature difference 2A7 = T hot - T cold is finite, the thermal 
rectification is obtained with AQ > in Fig. \T\ It means T{pi\,0)-i) > at these 7 's hence 
there is a greater heat flux in the direction from C reservoir to L reservoir than the heat flux in 
the inverse direction when the temperature bias is reversed. This indicates that by suffering the 
inelastic scattering introduced by the temperature probe, ballistic phonon has more probability to 
transport from C reservoir into L reservoir than transport in the inverse direction. To understand 
this, one must recall the fact that AQ is determined by the factors N(a>i, u) 2 , 7 ) and T((x>\, oj 2 ) as 
shown in Eq. ©. The phonons in reservoirs obey the Bose-Einstein statistics, as well as the fact 
that d{x[l + 2n(x)]}/dx — > when x — > as shown in Eq. ©, the phonons with the frequencies 
in the region < co < 37 /7 p /, (where the unit, k B T p h/h, is omitted) dominate N(o>i,cl>2, 7 ) when 
A7 is a small finite value and the phonon with the lower frequency is more dominant. Thus for 
Wt = 21 nm, at 7 = 0.37 p /,, the phonons with frequencies < co < 0.9 dominate N(to\, 002, 7 ). 
The phonons with frequencies < co < 0.25 is more dominant to N(co\,co2,Tq), but as shown 
in Fig. [2l t rc /t rl decreases much slowly in this region and it increases much faster in the region 



0.25 < co < 0.58. Hence by considering the both contributions from N(co\,co 2 , To) and T{to\, coa), 
the phonons with frequencies 0.25 < co < 0.58 determine AQ > because T(tO\, oi%) > in this 
region. The similar analysis can be applied to other TVs. One can expect that at very low T , 
AQ < because T(co\,co 2 ) < in the dominant region of frequency. It is confirmed by the inset 
ofFig.|Zlat7o = 0.17>. 

The thermal rectification can be studied by the efficiency which is defined as 77 = |AQ|/|<2-| x 
100%. The numerical results of the efficiency are shown in Figs. [3] and |4] for different TVs and 
different WV's. The efficiency can achieve about 5.5% as shown in the two figures. For the system 
with W T = 27 nm, the region of frequency in which t rc /t rl slowly decreases is the shortest 
as shown in Fig. [21 Thus at T = 0.3T ph , 0AT ph and 0.5T ph , this system achieves the highest 
efficiency. However, at T = 0.6T ph , the system with W T = 13 nm has the highest efficiency. This 
is because tr C /t rl increases in the region 0.35 < co < 0.7. 

In summary, we have studied the ballistic phonon heat flux flows in the asymmetric three- 
terminal systems and the ballistic thermal rectification is attained. The rectification is purely 
a quantum effect and is critically determined by the dependence of the ratio trc(co)/trl(co) on 
co. trc(co)/t rl (co) is determined by the inelastic scattering of phonons in the temperature probe. 
Therefore, one can attain the ballistic thermal rectification by introducing the inelastic scattering 
of phonons into the mesoscopic systems. 
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